Tight Noise Thresholds for Quantum Computation with Perfect Stabilizer Operations 
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We study how much noise can be tolerated by a universal gate set before it loses its quantum-computational 
power. Specifically we look at circuits with perfect stabilizer operations in addition to imperfect non-stabilizer 
gates. We prove that for all unitary single-qubit gates there exists a tight depolarizing noise threshold that 
determines whether the gate enables universal quantum computation or if the gate can be simulated by a mixture 
of Clifford gates. This exact threshold is determined by the Clifford polytope spanned by the 24 single-qubit 
Clifford gates. The result is in contrast to the situation wherein non-stabilizer qubit states are used; the thresholds 
in that case are not currently known to be tight. 
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Introduction. A way to study the resources needed for 
universal quantum computation (UQC) is to analyze the tran- 
sition from a system that can provide UQC to one that is clas- 
sically efficiently simulable. A particularly useful example of 
such a system is given by the stabilizer operations, which are 
made by a combination of: preparation of |0) states, unitary 
Clifford gates, measurements in the {|0), |1)} basis and clas- 
sical control determined by the measurement outcomes. The 
Gottesman-Knill theorem tells us that the such stabilizer op- 
erations can be efficiently simulated classically (see for exam- 
ple (J, Theorem 10.7]), while it also known that the addition 
of any other one-qubit gate outside the Clifford group will en- 
able the system to perform UQC. This fact provides us with 
a framework for testing tolerance to noise — one can examine 
how noisy this additional non-Clifford gate can be before it 
becomes classically simulable itself. If the non-Clifford oper- 
ation has become a probabilistic combination of Clifford gates 
due to the noise, then we know that we are unequivocally in 
the classical computational regime. The noise rate where the 
extra gate becomes simulable (where it enters the "Clifford 
polytope" flUl ) is an thus upper bound for fault tolerance. If 
the converse is true — i.e. if any operation outside the Clifford 
polytope enables UQC then the threshold is tight. In this arti- 
cle we show that for single-qubit gates undergoing depolariz- 
ing such a tight noise threshold does indeed apply. We will do 
so by proving that any depolarized gate that lies outside the 
Clifford polytope of single-qubit operations, in combination 
with noiseless stabilizing operations, allows for UQC. This 
result should be contrasted to the situation for non-stabilizer 
qubit states where the thresholds in that case are not currently 
known to be tight. 

We will consistently assume that Clifford gates can be im- 
plemented perfectly, motivated by the fact that these gates can 
be implemented fault-tolerantly by applying them transver- 



sally and to encoded states J2> d, B Htl- The fault-tolerant 
implementation of Clifford gates naturally carries with it a 
threshold of its own, independent of the kind we discuss in 
this paper. The current model is particulary relevant to the 
so-called Pfaffian quantum Hall state in topological quantum 
computation Il2l [l3|] , where the whole Clifford group (but 
only the Clifford group) can be implemented using braiding 
making these operations naturally fault-tolerant. The addi- 
tional resource required to perform UQC will likely be highly 
noisy, and so we can see the parallels with our model. 

We will begin by listing a couple of previously known re- 
sults in this area. Next, we will discuss the connection be- 
tween the geometry of the Clifford polytope and stabilizer 
measurements, and show that tightness of a magic-state distil- 
lation procedure for single-qubit states automatically ensures 
tight thresholds for non-Clifford gates undergoing any kind 
of unital noise. Finally we show that currently known magic- 
state distillation techniques are sufficient to prove tight thresh- 
olds for a non-Clifford gate undergoing depolarizing noise. 

Previously known results. The idea of using perfect stabi- 
lizer operations in conjunction with imperfect non-stabilizer 
states to perform UQC originates with Bravyi and Kitaev Jfl . 
The conditions on the ancillary qubits to enable UQC is that 
they are sufficiently close to being one of the 20 so-called 
magic states that lie on the surface of the Bloch sphere. The 
two classes of magic state (see Figure[TJ) are the \H) type and 
\T) type, where all \H) type states can be derived by apply- 
ing a Clifford operation to some canonical representative (and 
similarly for \T) states); 
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The routines used in JfJ] were unable to distill qubit states just 
outside the edges and faces of the octahedron of FigureQ] Re- 
ichardt |7j] subsequently suggested an improved routine that 
closed the gap in the \H) direction (along the edges of the 
octahedron). Virmani et al. U suggested using the convex 
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FIG. 1: Magic States and the Octahedron: Some of the single-qubit 
magic states: \H) type states are designated with black arrows, \T) 
type states with white arrows. The octahedron defined by |x| + \y\ + 
|z| < 1 depicts the single qubit states that can be created by stabilizer 
operations. Reichardt [7j] showed that distillation techniques work 
right up to the edges of the octahedron (i.e. tight in the \H) direction). 
Current distillation techniques are unable to distill states just outside 
the faces of the octahedron (i.e. not tight in the \ T) direction). 



hull of Clifford operations in order to find gates' robustness 
to various types of noise. In particular they considered gates 
that are diagonal in the computational basis. Plenio and Vir- 
mani 1 1 Qfl subsequently extended this idea by analyzing cases 
where noise was allowed to affect the stabilizer operations too. 
Buhrman et al. ifTltl used a similar idea (that noise causes non- 
Clifford gates to eventually become implementable via Clif- 
ford gates only) to find the non-Clifford gate that is most re- 
sistant to depolarizing noise — a 7t/8 rotation about the Z axis 
(or the same gate modulo some Clifford operation). Reichardt 
|8[] showed that this particular gate enabled UQC right up to 
its threshold noise rate (about 45%), as well as considering 
in detail the process of reducing multi-qubit states to single- 
qubit states using postselected stabilizer operations. Our cur- 
rent result here generalizes this tightness result to all possible 
single-qubit gates. 

Preliminaries and Notation. Let us parameterize an arbi- 
trary single-qubit SU(2) gate as follows 



£7(9,7,8) = 



e^cos(9) 
e- l5 sin(9) 



-e l8 sin(9) 
-^cos(e) 



(3) 



The representation of this rotation in SO(3) is denoted by 
R(Q,J, 8). Implementing a rotation R while suffering depolar- 
izing noise (with noise rate p), means that this noisy operation 
is represented by the rescaling M = (1 — p)R, a fact that we 
will need later. 

Often we will apply the unitary C/(8,y, 8) to one half of an 
entangled Bell pair, |<5) = -4=(|00> + 111)), yielding 



(/®t/)|4>)(4>|(/®t/) t 



(4) 



If we use the two-qubit Pauli operators as a basis for the 
density matrix p then we can find the 16 real coefficients 



Cjj = Tr (p (a,- ®Gj)) so that 

p = zEcyfa®^) *J e {i,x,y,z}. 



(5) 



Since we have applied a local unitary to a maximally entan- 
gled state, the coefficients (cix,C]Y,Ciz,cxi,Cyi,Czi) are al- 
ways zero. Comparing the 9 coefficients {cxx,cxy, ■ ■ ■ ,czz} 
one can see that these are the same as the entries of the SO(3) 
matrix R(Q, y, 8). More precisely, 

/ c X x 

R(B,y,S)=\ c XY -cyy c zy ) (6) 



CXX 
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czx 
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CXZ 


-c YZ 


czz 



where the c, 7 - are obviously also functions of (8,y, 8). 

If we represent the 24 single-qubit Clifford operations as 
SO(3) matrices, then they are simply signed permutation ma- 
trices with unit determinant (they are a matrix representation 
of the elements of the chiral octahedral symmetry group or, 
equivalently, the symmetry group S4). We label these opera- 
tions Cj and so the convex hull of the C, (the so-called Clifford 
polytope) is given by 




24 



with pi > and V pi = 1 



(7) 



Geometrically, the Clifford polytope is a closed polyhedron 
in M 9 that has 24 vertices (each vertex representing one of 
the Q). This polytope can also be defined by the bounding 
inequalities of its 120 facets. The concise description of these 
facets used by Buhrman et al. |l 1J is given by the set 



7 ={C,FC J \i 1 je{l,...,2^},Fe{A,A T ,B}} (8) 



where 




and B 



(9) 



At times we will have reason to refer to different subsets of 
the set of facets f so we use the obvious notation 



7 = TaU7 a tU7 B - 



(10) 



It is useful to note that all the facets derived from A comprise a 
single column with ±1 entries and zeros elsewhere, and simi- 
larly for the row facets derived from A T , hence \7a\ = \?a t \ = 
3 • 2 3 = 24. There are \7b\ =72 "B-type" facets, which can 
be constructed as follows: (i) Pick one position in the matrix 
e.g. Mj j and put ±1 there (9 x 2 = 18 choices), (ii) Fill the 
remaining entries not in row i or column j with ± 1 such that 
det(M) = -2 (4 choices). 

To determine whether or not an operation M is inside the 
Clifford polytope T we take the elementwise inner product 
(or Frobenius inner product) between M and the facets F £ 7 
of the polytope 



M ■ F = £ MijFij = Ti(M T F). 

i.i=i 



(11) 



Using the above notation, a 3 x 3 matrix M is inside the poly- 
tope i> if and only if for all F £ 7 we have M ■ F < I. 
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Interpreting the facets of the Clifford polytope. Our proof 
will involve applying some non-Clifford gate to one half of 
a Bell Pair (as in Eq. HI and then postselecting on the out- 
comes of various stabilizer measurements. This has the effect 
of taking our two-qubit state p to a single qubit state p' (times 
some stabilizer state that we do not care about), which we 
then distill using magic state distillation (see if for a more 
general discussion of these kinds of techniques). For exam- 
ple, performing a YX measurement on p and postselecting on 
a "+1" outcome (i.e. projecting with II = i (II + YX)) leads 
to a single-qubit state p' with a Bloch vector given by 



r(p') = ( 



Cll+Cyx ' Cu+Cyx 



(12) 



The form of this vector means that it lies in the YZ plane (see 
Figure Q}, where we know that distillation techniques work 
right up to the edge \y\ + \z\ = 1 of the octahedron. We can 
check if r is outside the octahedron by simply comparing the 
L l norm of r with 1. Rearranged, the condition ||?||i > 1 for 
being outside the octahedron is 



\CXZ - CZY | + | - (CXY + Czz) | > |c// + C Y X \ 



(13) 



Given the correspondence between the coefficients en and the 
elements of R (see Eq. [6]l we can rewrite the above condition 
(dropping the absolute value operators) as a facet inequality 



R ■ F > 1 where F 




(14) 




This facet is a legitimate "B-type" facet and a little thought 
shows that, had we applied the single qubit Pauli operations 
(as SO(3) rotations) X, Y or Z to r(p') above, we would arrive 
at three other "B-type" facets 



(15) 



respectively. Note that all four facet inequalities combined 
could be simplified to the form Eq. [13] above. We omit the 
details, but it should be easy to see that all 72 "fi-type" facets 
correspond to postselecting on some (weight two) Pauli oper- 
ator, and possibly performing a single qubit Pauli rotation on 
the resulting p'. 

It is somewhat more straightforward to see the geometri- 
cal interpretation of the "A ^ -type" facets. For example, the 
canonical A given in Eq. |9j merely returns the sum of the ele- 
ments of the Bloch vector r, arising from a rotation applied to 
the X eigenstate. 



RA 



3 
i=\ 



where r 




In general, an operation M having an inner product greater 
than one with some "A-type" facet simply means that M, ap- 
plied to some initial vector corresponding to a stabilizer state, 
brings that vector to a final position outside the octahedron. 



The preceding discussion shows us that if magic state dis- 
tillation was possible everywhere outside the octahedron, then 
every unital operation outside the Clifford polytope would be 
distillable — either straightforwardly or by using postselec- 
tion, depending on which facet it violated. Using current (not 
tight) distillation routines however, we would be unable to 
deal with some operations violating an "A-type" facet by a 
fairly small amount. In the next section we show that, for 
depolarizing noise, any noisy rotation violating an "A-type" 
facet also violates a "B-type" facet. Since "Z?-type" facets cor- 
respond to \H) state distillation, the results we obtain are tight. 

Tight threshold for depolarizing noise. The claim we shall 
prove is that, anytime a matrix M = (1 — p)R, representing a 
depolarized rotation, is outside some "A-type" facet then there 
exists a "B-type" facet that M also lies outside. In fact, we will 
prove the slightly stronger statement that for all R £ SO(3) 

VA G 7 A U 7 a t, 3B G 7b such that R ■ (B — A) > 0. (17) 

To simplify the proof we will repeatedly make use of the 
symmetries of the problem (see Eq. [gj. We will pick a canon- 
ical "A-type" facet and assume that this gives the largest inner 
product with R of all the F G 7a- If there was a larger inner 
product with some F G 7 a t, then we could just relabel R T as 
R. We can assume that the facet with ones in the first column 
(the A in Eq. |9]) gives the biggest inner product since 

Tr(R T {QACj)) = Tx(CjR T CiA) = Tr((C k RQ) T A) . (18) 

The proof will hinge on an entry of R, outside of the first 
column, being larger than the rest of the elements outside the 
first column. As such, let us define 12 matrices closely related 
to A and call them A' ; : 



a; = 




a'. 



12 




such that the index i largest inner product R ■ A\ tells us the 
sign and location of the largest magnitude element outside the 
first column. Once again, symmetry allows us to assume that 
Aj yields the largest inner product because the rest of the A\ 
can be derived from A\ via Clifford rotations 



{4} 




je {1,2,3} 
ke {1,2,3,4} 



(19) 

For a matrix R to be an SO(3) rotation there are constraints 
ngj on the signs of the elements Rj j i.e. there are 8 choices for the 
first column, 6 choices for the second column and 2 for the 
third. Given that A is the maximum facet for R, we have fixed 
the signs positively in the first column, reducing the number of 
types of rotation to 6 x 2 = 12. Since A\ gives the maximum 
inner product with R of all A' we have that R\ 2 < 0, which 
reduces the number of rotation types further to 3 x 2 = 6. It 
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FIG. 2: Proof of Eq.[23] For any pair of 2- vectors u and v with the 
same I? norm, the vector with greater L°° norm has smaller L 1 norm. 
A vector pointing towards the black dot has simultaneously minimal 
L°° norm and maximal L 1 norm. 

can be shown that R\2 having larger magnitude than the rest 
of the elements Ru(i £ {1,2,3},,/ £ {2,3}) restricts the type 
of rotation further to one the following four types 




This should not be surprising if one considers that 

#1,2 = -(#2,1*3,3 -#2,3#3,l) (20) 

because of the structure of SO (3) matrices, and the sign pat- 
terns listed above ensure |i?i,2| is as large as possible. 

We claim that the B £ fg of Eq.|9]will suffice to prove the 
desired inequality R ■ {B — A) > 0, which reads in matrix form 







i 
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>0. (21) 



Using the relevant entries of R we define a pair of 2-vectors u 
and v as 

U = (R l , U R L2 ) V = (/?2,3,#3,3) (22) 
then we can rewrite the above inequality Eq.l2Tlas 

||v||i-||5||i>0. (23) 

The L 2 normalization of all the rows and columns of the 
rotation matrix R means that u and v have the same L 2 norm. 
With reference to Figure |2j it should be clear that because u 
has an L°°-norm at least as big as that of v (because \R\j\ > 
|#2,3 1 , |#3,3 1), it holds that the L 1 -norm of v is automatically at 
least as large as the L 1 -norm of u, as desired. 

Summary. We showed that for any unitary 1-qubit gate 
undergoing depolarizing noise with rate p it holds that if its 
SO(3) representation M = (1 — p)R lies outside the Clifford 
polytope <E , then it must be the case that there a facet B £ J# 
such that M B > 1 . In turn, this means that if this noisy gate is 
applied to a Bell pair |4>) = -4=(|00) + |11)) and an appropri- 
ate stabilizer measurement is applied, then, conditionally on 
the outcome of the measurement, one obtains a state that can 
be transformed using Clifford gates into a single qubit state 
with |y| + \z\ > 1 in the Bloch ball representation. By the re- 
sult of Reichardt 0] such states enable stabilizing operations 
to perform universal quantum computation. 
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